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Abstract. We investigate the dynamics of a single deformable self-propelled particle 
which undergoes a spinning motion in a two-dimensional space. Equations of motion 
are derived from the symmetry argument for three kinds of variables. One is a 
vector which represents the velocity of the centre of mass. The second is a traceless 
symmetric tensor representing deformation. The third is an antisymmetric tensor for 
spinning degree of freedom. By numerical simulations, we have obtained variety of 
dynamical states due to interplay between the spinning motion and the deformation. 
The bifurcations of these dynamical states are analyzed by the simplified equations of 
motion. 
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1. Introduction 

Dynamics of self-propelled particles have attracted much attention recently from the 
view point of nonlinear science and non-equilibrium statistical physics. In non-biological 
systems, experiments have been conducted, for example, for oily droplets [U El [3j HJ |5] 
and Janus particles [El [7] which undergo chemical reactions with the molecules in 
the surrounding media. Theoretical studies have also been carried out by computer 
simulations both for the motion of an individual particle and for the dynamics of 
interacting particles [H [10] . Most of the studies assume that the particles are rigid 
without shape deformation. However, it is pointed out that, if a particle is soft, its shape 
is deformed when the migration velocity is increased as has been observed experimentally 
P] and analyzed theoretically [HI [12] . 

In biological systems such as living cells and micro-organisms, shape deformation 
plays a central role. In fact, migration of biological objects is induced by shape 
deformation and there are experimental investigations of the correlation of shape change 
and cell migration [T3[ [HI [151 EE E] • Theoretical studies for cell migration have been 
developed [TH1 HH1 [20l [21]. Computer simulations of self-propulsion driven by shape 
changes have also been available based on artificial deformable systems [221 123] . 

Recently we have introduced a model system for deformable self-propelled particles 
and have studied motion of individual particles and collective dynamics of interacting 
particles [211 [251 [261 [271 [28] . Dynamics of a single particle under external forces have also 
been investigated [29]. This set of model equations has been derived, by an interfacial 
approach, from an excitable reaction diffusion system both in two and three dimensions 

[301 EH- 

In this paper, we extend our model system to take into account a spinning motion. 
There are many experiments of spinning self-propulsion in both non-biological and 
biological microscopic systems, (i) One is an oily droplet to whicha small piece of solid 
soap is attached [32] . This causes a spinning motion as well as the ordinary translational 
motion of a droplet, (ii) The second example is a bacterium Listeria monocytogenes 
which causes locomotion together with spinning (helical) motion by polymerization of 
actin filaments [331 EH ESI [36] . Flagellated bacteria such as Escherichia coli also exhibits 
a spinning motion by rotating helical filaments [371 138] . In these biological examples, 
the left-handed and right-handed symmetry is broken because the filament has a specific 
rotating direction, (iii) The third is an anisotropic doublet composed of paramagnetic 
colloid particles undergoes a translation motion and rotation near a flat boundary under 
an oscillatory magnetic field [39]. The rotating direction follows the rotation of the 
magnetic field, (iv) The fourth is a discovery of spontaneous formation of spiral waves 
in a Dictyostelium cell, which is coupled with shape deformation, migration and rotation 
of the cell pEO]. 

In order to make our model as general as possible, we do not rely on any specific 
objects but derive the model equations only from the symmetry argument in terms 
of the antisymmetric tensor variable for spinning motion, the vector variable for the 
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translational motion and the traceless symmetric tensor for elongation of a circular 
particle. We keep the parity symmetry but allow a spontaneous symmetry braking in 
our model. Numerical simulations have been carried out in two dimensions to obtain 
variety of dynamical states caused by the interplay between the deformation and the 
spinning motion. We have also analyzed theoretically the bifurcations of these dynamical 
states. The case (ii) above is not considered in the present paper since the angular 
vector of filament rotation is parallel to the migration velocity and hence the dynamics 
is essentially three-dimensional. 

The organization of this paper is as follows; in the next section (Sec. [2]), we introduce 
the model equations. The results of numerical simulations are described in Sec. [3j The 
variety of the dynamical states are obtained as shown in the dynamical phase diagram. 
In Sec. EJ the theoretical analysis is carried out to derive the bifurcations of these 
dynamical states. Discussion for the results obtained is given in Sec. El 



2. Time-evolution equations 



We consider a self-propelled particle which changes its shape depending on the migration 
velocity. The degree of freedom of spinning motion around the centre of mass is also 
introduced. The basic dynamical variables are the velocity of the centre of mass Vi, 
the second-rank traceless symmetric tensor S^- for deformation and the antisymmetric 
tensor fly for spinning motion. By considering possible couplings among these variables 
and retaining some relevant nonlinear terms, the set of time-evolution equations is given 
by 
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and Oi, a 2 , bi, b 2 , b 3 , c\ and c 2 are the coupling constants and d is the dimensionality 
of space. Hereafter we will put d = 2. The repeated indices imply summation. The 
coefficient 7 takes either negative or positive values. Throughout the present paper, the 
parameters k and ( are assumed to be positive. The second rank traceless symmetric 
tensor Sij, which represents an elliptical deformation of a circular particle, is defined by 



St 
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where s > is the magnitude of deformation and the angle 9 represents the direction of 
elongation. The antisymmetric traceless tensor is defined as 



a 



o 



-id 



UJ 





(6) 



with u a dynamical variable. 
By putting (vi,v 2 ) = v (cos 
dv 
~dt 



,sin0), Eqs. (JTJ) - (JSJ) are written as 
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where ip = 8 — <p. Note that the second term on the right hand side of Eq. ([3]) vanishes 
in two dimensions. Because of the isotropy of space, only the relative angle if} enters 
into the set of equations as an independent variable. In fact, we have from Eqs. ([8]) and 

(HOD 



bi 



- — 1> 



a i 



s) sin 2ifj + {p2 — a 2 )uj- 



(12) 



dip 

~dt "2s' 2 

If the anti-symmetric tensor Q is not considered, the set of equations ([Q) and (|2J) 
has been studied both in two and three dimensions [2U [251 126]. Furthermore, those 
equations have been derived from a reaction-diffusion system [2DI EH]- See also Ref. 
[29] where dynamics of a self-propelled soft particle under external fields have been 
investigated. When the spinning degree of freedom is absent, i.e., a 2 = b 2 = b 3 = 0, 



Eqs. © - (Jin 



exhibit a bifurcation at 7 
kI + B 
2 B 



7 C where 



7c = 7T" 



with 



B 



"2k" 



(13) 



(14) 



When 7 < 7c, a straight motion of a particle is stable. However, this motion becomes 
unstable for 7 > 7 C and a circular motion (orbital revolution) appears. It should be 
noted that the coefficient in front of sin 2-0 in Eq. (fl2|) vanishes at 7 = j c . Another 
important parameter is the coefficient 61 in Eq. (J2J). When b\ is positive, a particle 
tends to elongate along the propagating direction whereas when it is negative, it deforms 
perpendicularly to the translational velocity [24J. Actually we have found that b\ < 
for the excitable reaction-diffusion system [301131]. Throughout this paper, we will put 
b\ < and a\ < so that the constant B defined by Eq. (fl4|) is positive provided that 
the relaxation rate of deformation k is positive. 
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Figure 1. (color online) Dynamical phase diagram obtained by solving Eqs. (0 - 
(fTTjl numerically. The meaning of the symbols is given in the text. The bifurcation 
boundary from the motionless state to the spinning state without migration is indicated 
by the broken line which is obtained from Eq. (|25j) . The thin solid line obtained from 
Eq. (|29|) represents the boundary between the v = and the u ^ region. The dotted 
line given by Eq. (j39| is the threshold at which the state u> = becomes unstable. The 
parameters are k = 0.50, a\ = — 1, a-i = 0.75, b\ = —0.5, &2 = 1, C — 2, and ci = 2. 
The bifurcation threshold 7 C defined by (|13p is given by 7 C = 0.75. 



Before closing this section, we summarize the meaning of other terms in the time- 
evolution equations ([[]), © and (EJ). It is evident from Eq. (JSJ) that the term with the 
coefficient a 2 tends to curve the trajectory of migration. We choose 02 > so that a 
particle spinning to the counter-clockwise direction (viewed from the top) rotates to the 
left (viewed from behind). This is similar to the Magnus force [H] which is proportional 
to the density of the surrounding fluid and the cross section of a rigid cylinder and the 
volume of a rigid sphere. The 5 2 -term in Eq. fl2]) is the same as the convective term 
for the orientational tensor in liquid crystals under rotational flow [42J so that we put 
b 2 = 1. In other words, the 6 2 -term is not dissipative whereas all other terms on the right 
hand side of Eq. fl2]) are dissipative. The last term in Eq. (J2J) with b 3 > represents 
that spinning motion enhances deformation whereas the C2-term in Eq. ([3]) with c 2 > 
has an effect such that an elongated particle prevents from spinning as can be seen 
from Eq. ( fill) . The ci-term drops out in two dimensions. The coefficient ( in Eqs. (Hj) 
is chosen to be positive. This means that the particle has an internal mechanism to 
keep spinning motion. Furthermore, the potential G takes a form that both clockwise 
and counter-clockwise motions are equally possible. It is noted that the coefficients in 
the v 3 term in Eqs. ([1]) and the trf2 4 term in (j2J) can be eliminated, without loss of 
generality, by absorbing those in the definitions of V{ and u. Quite recently, we have 
derived Eq. ([1]) with a\ = a 2 = in fluids taking into account the Marangoni effect and 
the hydrodynamic interaction [43]. However, such a study for shape deformation and 
spinning motion has not been available at present. 
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3. Numerical Results 

In the preceding section, we have fixed the sign of the coefficients in the time-evolution 
equations ([I]), ([2]) and (j3J) with (j4]) by considering the mechanisms of each term. 
However, it is difficult to determine the magnitude of those coefficients within the present 
phenomenological approach. We are concerned with the effects of spinning on the self- 
propelled motion. Therefore, we choose the coefficient 7 in Eq. (JT|) as an important basic 
parameter since three different kinds of solutions, motionless state, straight motion and 
circular motion, are realized by changing 7 in the absence of spinning. The other aspect 
to be considered is the coupling between spinning and deformation. In Eq. ([2]) for the 
deformation tensor, the coefficient b 2 has been uniquely put to be b 2 — 1- Therefore the 
remaining coefficient 63 of the nonlinear coupling between S and Q is chosen as another 
basic parameter. Other remaining coefficients are set to be of the order of unity as 
k = 0.5, a\ = —1, a 2 = 0.75, b x = —0.5, b 2 — 1, C — 2, and c 2 = 2. 

Figure Q] show the dynamical phase diagram obtained from Eqs. (J7J) - ( TTTi) . In 
numerical simulations, we have employed the fourth Runge-Kutta method with the 
time increment 5t = 10~ 4 . The different symbols correspond to different dynamical 
states as explained in detail below. 

The state indicated by the crosses in Fig. [1] is the motionless state. In this state, 
v = and s = and the value of u is finite. That is, the particle with a circular 
shape is spinning and its centre of mass does not move. In the region indicated by the 
filled circles, the variables s and u take finite constant values while the magnitude of 
the velocity is zero, v — 0. The angle of the elongated direction 9 varies monotonically 
in time. Therefore, the particle is elliptically deformed while the centre of mass does 
not move and undergoes clockwise or counter-clockwise rotations depending on the 
initial condition. We call this motion as a spinning motion throughout this article. 
The white square symbols and the plus symbols mean two different circular motions, 
which we call the revolution I and revolution II states, respectively. In both of these 
states, all of the variable u, s, and v take finite constant values and the angles <f> and 9 
varies monotonically in time, keeping their difference if) = 9 — <f> constant in time. The 
trajectories of the revolution I and II motions in real space are displayed in Figs. Efa) and 
(b), where some snapshots of the particle which undergoes counter-clockwise rotation 
are shown. In both states, clockwise rotation can also appear depending on the initial 
condition. It is noted in Fig. [T]that there is a region where the revolution I and revolution 
II states coexist. The properties of the revolution I and II states will be discussed in 
detail in Sec. HI 

There is another orbital revolution which appears for large values of 7 in the region 
indicated by the diamonds in Fig. [TJ We call this state as a circular state where u = 
while the values of v and s are finite constants and the angle of the propagating direction 
and the elongation direction 9 varies monotonically with a fixed finite difference 
if) = 9 — <f). The trajectory and the particle shape of a counter-clockwise circular motion 
in real space are displayed in Fig. EJ^c). Since no spinning motion occurs, this is nothing 
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Figure 2. (color online) Rotation and orbital revolution of a particle in the counter- 
clockwise direction. Trajectory in real space is displayed for (a) revolution I motion for 
7 = 1 and &3 = 0.5, (b) revolution II motion for 7 = 1 and 63 = 0.9, and (c) circular 
motion for 7 = 3 and b 3 — 0.5. The particle size in (a) - (c) is reduced by the factor 
of 1/4 for the sake of clarity, (d) Radius of the circular orbit in real space, (e) the 
relative angle ip, and (f) the spinning variable w as a function of 7 for 63 = 0.5. 



but the circular motion obtained in Ref. |24J. 

The quantitative differences of these three motions are summarized in Figs. EJ^d), 
(e) and (f). In Fig. Efd), the radius of the circular orbit in real space of the revolution 

I motion, the revolution II motion, and the circular motion is plotted as a function of 
7 by the squares, the circles, and the triangles respectively. In drawing a particle, the 
radius of an undeformed circular particle is set to be r = 2 throughout the present 
paper. The radius of the revolution I motion in a larger scale is also plotted by the open 
squares. Here, note that the radius of the revolution I motion becomes extremely large 
around 7 = 2.33. This will be analyzed theoretically in section HJ In Figs. and (f), 
the values of the variables ip and uj are plotted for u; > 0, respectively. Note that ip 
for the revolution I motion (the stars) is larger than ir/2 whereas i/j for the revolution 

II (the pluses) and circular (the crosses) motions is less than tt/2, which is consistent 
with the particle configuration in Fig. EJa). The values of if) of the revolution II state 
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Figure 3. (color online) (a) Trajectory and (b) return map of the quasi-periodic I 
motion for 7 = 1 and &3 = 0.8 in the real space. The trajectory for a shorter time 
interval is indicated by the thick solid line in (a). The return map is obtained at the 
Poincare section indicated by the horizontal line in (a). 



(a) 



2 

^ 





4 




-v1 


3V 









-2 2 

x 




Figure 4. (color online) (a) Trajectory and (b) return map of the quasi-periodic I 
motion for 7 = 1 and &3 = 0.9 in the real space. The trajectory for a shorter time 
interval is shown by the thick solid line in (a) together with the chronological order 
1, 2, 3, and 4. The return map is obtained at the Poincare section indicated by the 
horizontal line in (a). 



smoothly connect with those of the circular state at around 7 = 2.3. There is no such 
a connection with that of the revolution I state. Similarly, in Fig. [2]^f), the values of oj 
of the revolution II state (the triangles) continuously connect with those of the circular 
motion (the diamonds) at around 7 = 2.3, but no such a connection occurs with the 
revolution I state (the circles). 

In the region of the white circles around 7=1 and 63 = 0.75 in Fig. [TJ a quasi- 
periodic I state appears, whose trajectory of a counter-clockwise rotation is displayed 
in Fig. [3]^a). In Fig. [3]^b), we show the return map of this motion, which is obtained at 
the Poincare section indicated by the horizontal line in Fig. |3]Ja). In this state, all of 
the variables v, <f>, s, 6, and u, as well as the difference if) = 9 — <f> are time-dependent. 
There is another quasi-periodic motion called a quasi-periodic II state in the region of 
the triangles in Fig. [TJ One example of the trajectory in real space is shown in Fig. H^a) 
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Figure 5. (color online) Attractors in the s-tj) plane for the revolution I and II states, 
the quasi-periodic I and II states, the period-doubling state, and the chaotic state for 
positive values of ui and for 7 = 1 and (a) b 3 = 0.3, (b) b 3 = 0.5, (c) 63 = 0.8, (d) 
b 3 = 0.82, (e) b 3 = 0.844, and (f) 6 3 = 0.9. The broken line is the nullcline for ds/dt = 
and the dotted line is the nullcline for dtp/dt = obtained from the reduced equations 
(|3"Tjl and (|32"f . The double circle (white circle) is the stable (unstable) fixed point 
corresponding to the revolution I state. The black square indicates the saddle point 
and the filled circle indicates the stable fixed point corresponding to the revolution II 
state. The limit cycles in (c) and (f) correspond to the quasi-periodic I and II states 
respectively. The attractors in (d) and (e) correspond to the period-doubling state and 
the chaotic state respectively. 



for 7 = 1 and 6 3 = 0.9, where a particle rotates in the counter-clockwise direction. 
Some parts of the trajectory are highlighted by the thick solid line together with the 
number indicating the chronological order. In this state, all of the variables u and v 
and s are non-zero, and both of <f> and as well as their difference tp = — <fi are time- 
dependent. In Fig. 11(b), we show the return map of the quasi-periodic II motion, which 
is obtained at the Poincare section at the horizontal line in Fig. [4(a). Note that both 
the quasi-periodic I and II states coexist with the revolution II state (plus symbols) in 

Fig.m 

In order to elucidate further each motion described above, we consider the dynamics 
on the s-ip plane. Figure [5] shows the attractors for 7 = 1 and for different values of 63. 
In Figs. [5] (a) for 63 = 0.3 and (b) for 63 = 0.5, the stable fixed point indicated by the 
double circle is the solution of the revolution I state. In Fig. [5(b) another stable fixed 
point of the revolution II state indicated by the black circle appears as a saddle-node 
bifurcation. In Fig. [5(c) for 63 = 0.8, the fixed point corresponding to the revolution I 
state becomes unstable via a Hopf bifurcation and a limit cycle oscillation appears which 
represents the quasi-periodic I state. Increasing 63, the simple limit cycle of the quasi- 
periodic I state becomes unstable and a period-doubling occurs as shown in Fig. [5(d) 
for 63 = 0.82. By increasing 63 further, a chaotic state appears as shown in Fig. [5(e) 
for 63 = 0.844. When 63 is extremely large, it turns out that ip varies from to 7r as 
shown in Fig. [5(f) for 6 3 = 0.9, which corresponds to the quasi-periodic II state. Figure 
[5] shows the case of positive values of u. The attractors and the nullclines for negative 
values of uj can be obtained by replacing ip by n — ip. 
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Figure 6. (color online) Trajectory in the s-cj) plane for (a) the quasi-periodic I state, 
(b) the period-doubling state, (c) the chaotic state, and (d) the quasi-periodic II state. 
The thicker (dotted) lines indicate the trajectories for a shorter (longer) time interval. 
The parameters are the same as those in Fig. (5][c) - (f ) respectively. 
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Figure 7. (color online) Trajectory in the real space of (a) period-doubling state and 
(b) chaotic state. The trajectory for a shorter time interval is shown by the thick solid 
line. The parameters are chosen as 7 = 1 and (a) b 3 — 0.82 and (b) 63 = 0.844. 



Now a question arises. Why does a limit cycle oscillation in the s-ip plane cause the 
quasi-periodic I motion in the real space? Similar problems exist for the period-doubling 
too. The origin can be traced back to the isotropy of space. In fact, the dynamics of 
the system is governed by the set of Eqs. (jTJ), ©, jUJ) and (TT21) whose solutions behave 
as shown in Fig. [5j However, when we plot the trajectory of the particle in the real 
space, we have to solve Eq. (J8]) for the angle of the velocity. The magnitude s of 
deformation is shown in Fig. [6] as a function of <ft for the parameters corresponding to 
those of Fig. [5(c) - (f). It is evident from Figs. [5(a) and [6(a) that the trajectory on 
the s — ip — space composes a torus and hence a quasi-periodic I motion can appear. 
Other complex dynamics of the period doubling, chaotic motion and quasi-periodic II 
motion can also be understood from Figs. [6(b), (c) and (d), respectively. 

As explained above, the limit cycle corresponding to the quasi-periodic I state 
becomes unstable by increasing the values of 63 and the period-doubling bifurcation 
occurs as shown in Fig. [5(d). The corresponding trajectory in the real space is displayed 
in Fig. [7](a). Further increase of 63 leads to a chaotic motion as shown in Fig. [7(b). This 
chaotic state eventually disappears for larger values of 63 and, in turn, the quasi-periodic 
II state appears. This series of the dynamical transitions from the quasi-periodic I state 
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Figure 8. Bifurcation diagram near and in the chaotic regime obtained from Eqs. ([7]) 
- (fTTj) with 7=1 imposing the conditions cos 2-0 = 0.7, sin 2-0 < 0, dip/dt > 0, and 
uj > 0. 



for 63 = 0.8 to the quasi-periodic II state for 63 = 0.9 are not shown in Fig. [T]due to 
the size of the grid. However we show, in Fig. |8l the bifurcation diagram in the vicinity 
of the chaotic regime is displayed, which has been obtained numerically from Eqs. (|7|) - 
( TTTT) with 7 = 1 under the conditions cos20> = 0.7, sin 2ip < 0, dip/dt > 0, and uj > 0. 
This is a typical route to chaos via period-doubling bifurcation well known in a logistic 
map [H]. Because of the choice dip/dt > 0, the quasi-periodic II state with oj > 0, 
where ip decreases monotonically, does not appear in Fig. |HJ 

4. Analysis of bifurcations 

In this section, we analyze the bifurcations found numerically in Section |3j However, 
direct analysis of set of Eqs. - ( 1T2|) for the four independent variables seems 
impossible. Therefore we have to introduce some approximations or simplification. Since 
our main concern is the coupling between spinning motion and deformation, we retain 
the deformation degrees of freedom 5* and eliminate the variables v and u by putting 
dv/dt = d/u/dt = in Eqs. ((7|) and (TIT]) , respectively. The stability and bifurcations 
of the solutions to the set of Eqs. ((9]) and ( 1T2|) are investigated. It is emphasized that 
this approach is justified when the particle is sufficiently soft and near the stability 
threshold where a rectilinear motion becomes unstable in the absence of spinning, that 
is, for k <C 1 and 7 ~ 7 C where 7 C is given by Eq. ffl3|) . 

Numerical simulations in the preceding section have been carried out in the 
condition that all the parameters are of the order of unity. Therefore the comparison 
with the predictions obtained by the reduced set of Eqs. ([9]) and ([12]) is expected to 
be qualitative. Nevertheless, it will be found that some of the phase boundaries in 
the dynamical phase diagram in Fig. [1] agree semi- quantitatively with the theoretical 
results. 
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By putting dv/dt = in Eq. (j7|), the velocity v can be expressed as 
when T < 



where 



T 1 / 2 whenr>0 ' (15) 



r = 7 — ^-s cos 2ip. (16) 



Equation ffTTj) with du/dt = gives us three stable solutions: 




i " forC-c 2 s 2 <0 

— forC-c 2 s^>0 ' (17) 



where 



a( 8 ) = (c-c2^r ■ (is) 

First, we consider the motionless state v — 0. In this case, the angle of the 
propagating direction <p has no meaning. Therefore, Eqs. and (fTUj) become 

— = — KS + D3SUJ , (19) 

§ = M . (20) 

Since uj is given by Eq. (iT7|) . Eq. (|T9"|) is closed for s. Therefore, it is sufficient to analyze 
the stability of the fixed point of Eq. ( TT9"j) . If oj = 0, Eq. (IT91 becomes, 

ds , , 

— = -ks. 21 

Since /t > 0, we have the solution s = 0. However, this gives rise to ( — c 2 s 2 = ( > 0, 
which contradicts the condition ( — c 2 s 2 < for u> = in Eq. (Il7p . On the other hand, 
if uj = ±u), Eq. f fl9|) becomes 

J = -& 3 c 2 (, 2 - if) 8, (22) 

where we have defined 

K = (23) 
c 2 

« = -. 24 

h 

We obtain a stable solution s = for K < 0, which loses its stability at if = and 
a pair of other stable solutions s = ±\[~K appears via a pitchfork bifurcation. The 
stability condition ( — c 2 s 2 > of uj = ±d) in Eq. ( !T5|) is satisfied for both of the stable 
solutions. The rotating frequency u is given from Eq. ( llSp by iy^ for s = and 
for s = ±\/if , respectively. 

We identify these stable fixed points with the states with v — obtained from 
the original equations (J7|) - ( II First, note that, when s — 0, the circular particle is 
not deformed and hence the angle of the longitudinal axis of the deformation 9 loses 
its meaning. Therefore, the solution (s, uj) = (0, ±a/C) corresponds to the motionless 
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state. In the same way, the pair of the solutions (VF, ±a/£) corresponds to the counter- 
clockwise and clockwise spinning motion. The pitchfork bifurcation mentioned above is 
a bifurcation from the motionless state to the spinning state. The bifurcation threshold 
6 3 = 63 is given by the condition K = as 

b; = - c (25) 

The bifurcation threshold (125]) is indicated by the thick broken line in Fig. [1] in a good 
agreement with the numerical results obtained from Eqs. ([7]) - f fTT]) . 

The condition T < for v = in Eq. ( fl5l) should be satisfied for both the 
motionless state and the spinning state without migration. In the motionless state, 
which is represented by the stable fixed point (s,u) = (0, ±a/C), the condition becomes 
T = 7 < 0. On the other hand, the derivation of the condition for the spinning state 
in terms of the original parameters is more involved. As we have shown above, the 
variables are given by (s,u) = (VK, in the spinning state. By substituting these 

variables into Eq. ( fl6|) . the condition in Eq. ( fl~5|) is written as 

r = 7 - yfsT 1/2 cos2^ < 0. (26) 

Although the variable ip loses its meaning when v = 0, we may define it around the 
bifurcation from v — to v ^ 0. In fact, from Eq. ffT2~]) with v = 0, we obtain 

^ = ^K^sin2^ + X Ub 2 -a 2 )R^, (27) 
where Xu is the sign of u. The steady solution is given by 

^ = - 2x ^-^ (I) 1 ' 2 . (28) 



ai \K / 

Here, the stability condition of the solution ( 128]) is given from Eq. ( 1271) by a± cos 2ip < 0. 
Then, by using Eq. ( 1281) . the condition (T2"6"]) for v = becomes 7 < 7^, where 

7 t=|(^-i^^y /2 . (29) 

This expression is valid as long as Eq. ( 128]) satisfies | sin 2ip | < 1 which imposes a 
constraint for 63 as 63 > 63, where 

6 t = ^{ 1 + 4e 2 fep)* j (3()) 

On the other hand, for 63 < 63 < 63, Eq. f[2T]) has no steady solution, and hence ip varies 
from to 71. Therefore, on the average of cos 2ip, the condition for v = 0, given by 
Eq. ( 126]) . becomes 7 < 0. Consequently, the stability threshold of v — in Eq. ( fl~5]) is 
given by 7 = for < 63 < b\ and 7 = 7^ for 63 > 63 where 63 = 0.375. The bifurcation 
threshold from v = to a finite v obtained in this way is shown by the thin solid line 
in Fig. [T] which is consistent with the numerical results of Eqs. ([7]) - (TTT]) . 
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Next, we derive the stability condition for the migrating states v 7^ 0. The velocity 
is given from Eq. (j!5p by v — VT with T given by Eq. ffTB"]) . Then, from Eqs. ([8]) - ([TO]) , 
the time-evolution equations become 

_ = - K h + Bcos 2 2ij))s + 6i7cos2^ + o 3 sw 2 , (31) 
at 

-Bcos2ip s \ sm2ip + (02 — a,2)uj, (32) 



dt V 2s 2 2 



where w and B are given by Eq. (1171) and Eq. ( lT4j) . respectively. When oj = ±u), there 
is a constraint £ — C2S 2 > as well as s > and T > 0. The linear stability matrix of 
Eqs. flSU) and ([32]) is defined by 



L(s,i/>) 



(33) 



d s (ds/dt) djp(dsfdt) 
d s (dip/dt) d^dip/dt) 

where d s means the partial derivative with respect to s. The method we employ to obtain 
the stability of the fixed points is as follows; first, we evaluate numerically the fixed 
point of Eqs. ( 13~TT) and ( 132]) with oj = ±d) under the constraints T > and £ — c^s 1 > 0. 
Then, by using the determinant and the trace of the linear stability matrix ( 133]) . we can 
examine the stability of the fixed point. Note that the reduced equations ( |3TT) and ( 1321) 
are invariant under the simultaneous transition oj — > —oj and ip — > tt — ip. Therefore, we 
will consider only the case of uj = +Cj hereafter. 

The nullclines and the fixed points obtained from Eqs. ( 13T1) and ( 1321) with a; = a) are 
plotted on the s-ip plane as shown in Fig. |5] where the attractors obtained numerically 
from Eqs. (J7J) - ( fill are also superposed. The s- and -^-nullclines are displayed by the 
broken line and the dotted line respectively. The black squares and the white circles 
indicate the saddle points and the unstable fixed points of Eqs. ( 13~T1) and ( 1321) . The double 
circles and the black circles are the attractors of the revolution I and II states respectively 
obtained from Eqs. ([?]) - ( ITTj) . which agree, within numerical accuracy, with the stable 
fixed points obtained from Eqs. ( 13~1T) and ( 1321) . The limit cycle orbit corresponding to 
the quasi-periodic I state obtained from Eqs. ([7]) - ([TT]) exists around the unstable fixed 
point. The stability analysis of the fixed points of the reduced equations, Eqs. ( 13T1) and 
( 132]) with uj = Co, show a good agreement with the bifurcation behavior obtained from 
the numerical simulations of the original equations ([7]) - ( fTT]) . 

The stable and unstable blanches of the steady solution s for oj 7^ are displayed 
in Fig. Ufa) for 7 = 1 and by varying the value of 63. There are at most three fixed 
points; One is plotted by the thick solid line, which corresponds to the stable revolution 
I state. This fixed point loses its stability at around 63 = 0.66 by a Hopf bifurcation 
and becomes an unstable fixed point, displayed by the thick dotted-broken line. At 
around 63 = 0.48, a pair of another stable fixed point and a saddle point appears by a 
saddle node bifurcation, which are displayed by the thin solid line and the broken line, 
respectively. This stable fixed point corresponds to the revolution II state. 

Finally, we investigate the dynamics around the region oj = 0. The fixed point of 
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Figure 9. (color online) Bifurcation diagram varying (a) 63 for 7 = 1 and (b) 7 for 
63 = 0.5. The stationary value of s is obtained from the reduced equations (I5T1) and 
(|32|) with uj given by Eq. (fT7| . The thick and thin solid lines indicate revolution I and 
II states respectively. The dotted broken line and the dotted line indicate the unstable 
solutions. In Fig. (b), the stable solution with uj = 0, which corresponds to the circular 
state, is also shown (the line indicated by "c"). 



Eqs. ( l3lT) and ( )32l) with u = is given by 

617 cos 2ip 



k (1 + B cos 2 2^0 



(34) 



and 



tan2^ = for 7 < 7c, (35) 



C082 ^= 2g( 7 - k/2) fOT7 - 7c ' (36) 

where 7 C has been defined by Eq. (TT5j) . From Eq. (1511) . the requirement of the positivity 
of v in Eq. (JT3]1 gives us 

_ E7 cos 2 2?/; = 7 /o 7 n 

7 l + 5cos 2 2^ 1 + 5 cos 2 2iP ' 1 j 

Since B is positive from the definition ( JTlj) . this requirement is satisfied as long as 7 > 0. 
These results are consistent with those in Ref. [24] . where the spinning motion was not 
considered. The condition ( — c 2 s 2 < must be satisfied for 00 = in Eq. (1T7|) . From 
Eqs. ( )34|) and ( )36|) . this condition is written as 

_ c 2 (fe l7 ) 2 cos 2 2^ _ ^(7-^/2) 

^ /t 2 (l + 5cos 2 2V>) 2 ^ 2k£ " 1 ' 

This leads to 7 > 7* for c 2 > where the bifurcation threshold 7* is given by 
k 2kB( k CLi( 
2 + 6 2 c 2 ~2 + W 

It is noted from Eqs. flU and ([39]) that 7 C = 0.75 and 7* = 2.25. 

These results are summarized in Fig. [9](b), where the fixed point values of s and 
its stability are plotted for oj = together with those for oj = uj by changing the values 
of 7 for the fixed value of 63 = 0.5. For large values of 7, there is only one stable fixed 
point with u> = shown by the thin solid line in Fig. [9](b), which corresponds to the 



7* = 0+ ^= 0+ 77^- ( 39 ) 
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Figure 10. (color online) Angular velocity d<fi/dt, given by Eq. <JSj> , of the revolution I 
state with ui > obtained numerically from Eqs. (JT)) - (|TT|) . Between 7 = 2.33 and 2.34, 
the sign of d<f>/dt changes from positive to negative, where the radius of the circular 
trajectory in real space diverges as shown in Fig.[2jd). 




circular state. Around 7 = 2.427, another stable fixed point and an unstable saddle 
point with finite u appear by a saddle node bifurcation, which are shown by the thick 
solid line (re-I) and by the thick broken line, respectively. By decreasing the value of 
7, the stable fixed point, which corresponds to the revolution I state, becomes unstable 
around 7 = 0.8. On the other hand, the line of the saddle point crosses with the line 
of the stable fixed point of the circular state at around 7 = 7*. Then, the circular 
state loses its stability and the saddle point become stable by a transcritical bifurcation. 
This stable fixed point indicated by the solid line (re-II) corresponds to the revolution 
II state. 

It is mentioned that the anomalous behavior of the radius in the revolution I motion 
at 7 ~ 2.33 shown in Fig. [2]^d) is reproduced by the reduced set of equations (I3~TI) and 
fl32l) with Eq. (IE]). In Fig. [TUl we have plotted d(j)/dt as a function of 7 substituting 
the stationary values of the revolution I state. It is found that the angular velocity 
changes its sign around 7 ~ 2.33 from positive to negative by increasing 7. This causes 
an anomalously large radius of the circular trajectory in the real space and the change 
of rotating direction. 

In summary, the reduced equations ([311 and (1321) for s and ip, and the time-evolution 
equation ([8]) for <fi together with v = vT 7 and u = ±u reproduce all of the revolution I 
and II states, the circular state and the quasi-periodic I and II states. However, they 
do not reproduce the periodic-doubling state and the chaotic state. 

5. Discussion 

In this paper, we have investigated in two dimensions the dynamics of a deformable 
self-propelled particle having a spinning degree of freedom. By numerical simulations 
of Eqs. dl]) - ([3]), we have obtained a dynamical phase diagram as shown in Fig. [I] for 
k = 0.5, which displays a rich variety of dynamical states. Apart from the motionless 
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state, there are the spinning state, three types of revolution states and two quasi-periodic 
states. The trajectories in the real space of the latter five states are displayed in Fig. |2l 
Fig. |3]^a) and Fig. H^a). The period-doubling and chaotic motions, as shown in Fig. [7J 
are also obtained between the quasi-periodic I and II states. 

Theoretical analysis has been developed for the bifurcations between these 
dynamical states based on the two- variable equations f|T9|) and (1201 for v = 0, and 
Eqs. (131)) and fl32|) for finite values of v. This simplified set of equations succeeds 
in reproducing all the bifurcations except for the chaotic behavior between the quasi- 
periodic I and II states. Although not described here, we have verified that all of the 
dynamics can be reproduced if we employ the three- variable system in terms of s, ip, and 
bj. However, another set of three variables s, xp, and v is unable to realize the chaotic 
state. 

We have found numerically an anomalous increase of the radius in the revolution 
I state as shown in Fig. |2]^d). This anomaly originates from the existence of zero in 
d<p/dt as a function of 7 in the revolution I state as shown in Fig. [TQJ This is obtained 
numerically from the original equations ([7]) - (ITT)) for 63 = 0.5 with u > 0. The angular 
velocity d<p/dt is positive for 7 < 2.33 while is becomes negative for 7 > 2.34. Therefore, 
a particle in the revolution I state with u > undergoes a counter-clockwise orbital 
rotation for 7 < 2.33, whereas it rotates to the clockwise direction for 7 > 2.34. This 
implies that if a noise term is added in Eq. (JH)), a particle changes randomly the rotation 
direction in the vicinity of the anomalous point. 

Our representation of the basic set of equations (DO) - (BJ is independent of the 
dimensionality. Therefore, the present study of spinning motion of a soft particle can 
be readily extended to three dimensions. It is of particular interest to investigate the 
case that the spinning axis is parallel to the migration velocity. This corresponds to the 
motion of flagellated bacteria [371 EH] • We will return to this problem somewhere in the 
near future. 
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